function pdop = pdop(G)
    %% conpute the PDOP given the joint measurement matrix G
    pdop = trace(inv(G'*G));
end